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In a recent publication [PRE 86, 04012 (2012)], Santos has presented a self-consistency condition 
that can be used to limit the possible forms of Fundamental Measure Theory. Here, the direct 
correlation function resulting from the Santos functional is derived and it is found to diverge for all 
' densities. 
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In a recent contribution, Santos introduced a novel argument aimed at eliminating a source of ambiguity in the 
derivation of the Fundamental Measure Theory (FMT) approach to Density Functional Theory for hard spheres [![. 
The result is a new ansatz for improvement of FMT beyond the basic Rosenfeld functional [![• The proposal is quite 
interesting as the most accurate density functionals currently in use (such as the "White Bear" functional^ 0]) 
are of exactly this type: heuristic improvements beyond functionals based on Rosenfled's original reasoning together 
with the additional requirement that the forms reproduce known, exact results in low-dimensional systems. The 
introduction of a new element that eliminates some of the arbitrarity of these extensions is therefore welcome. The 
proposal of Santos is based on an exact scaling relation of the type successfully exploited by him and co-workers in 
the development of highly accurate approximations for the free energy and pair-distribution function of mixtures of 
hard spheres 0, Q. The purpose of this Comment is to examine one consequence of the proposed ansatz, namely the 
implied direct correlation function (DCF). 

The direct correlation function is a fundamental element in DFT as it provides a connection between model free 
energy functionals and liquid-state properties, for which much is known 7}. Given a (grand-canonical) free energy 
functional, D, [p] — Fid [p] + F ex [p] — pp, where p (r) , is the ensemble-averaged local density, Fa is the ideal gas 
contribution, which is not relevant here, p is the chemical potential and F ex [p] is the excess term, the (two-body) 
direct correlation function is given by taking two functional derivatives with respect to the density, 



6 2 pF ex [p] 

c 2 (ri, r 2 ) = - 1 J , 1) 

Sp (rx) Sp (r 2 ) 

where j3 = l/ksT , ks is Boltzmann's constant and T is the temperature 0, H| ■ This relation between the free energy 
functional and the DCF has always provided an important connection between free energy models and liquid-state 
properties: for example, one of the first indications of the utility of the White Bear functional was its improvement 
in the predicted DCF of hard spheres^. 

In DFT, the only unknown is the excess term and FMT is based on an ansatz of the form 



(5F ex = J $(n(r;M))dr (2) 
where the weighted densities have the generic expressions 

ni(r;[p]) = / Wi (r - r') p (r') dr' (3) 



Different models involve different collections of density- independent weight functions, Wi, and of different forms for 
the function $(n). The proposal of Santos makes use of the weight functions as were introduced by Rosenfeld 
(w s (r^) = 8 (§ — ria), w v (1*12) = O (| — ^12), w Vi (ri2) = f\2 t iS (% — r^) where a is the hard-sphere diameter) and 
the $ function of the Rosenfeld, = s<I>i (77) + $2 (v) (s 2 _ v 2 ) + ^3 ('?) s (s 2 _ 3w 2 ) , where 77 (r; [p]) — n v (r; [p]) is 
the weighted density formed from the weight function w v (1*12), etc. The other terms are 

$1 (V) = In (1 - V) , *2 (V) = ^(i^j, (4) 

The Rosenfeld functional reproduces the Percus-Yevik equation of state for a uniform fluid (p (r) = p with p being a 
position- independent constant) 0, 0] (see the Appendix for details). The idea of the extension discussed by Santos is 
that one would like to use knowledge of more accurate equations of state than Percus-Yevik to construct potentially 
more accurate approximations for <E>. Note that for a uniform system, Eq.(2) shows that the excess free energy of 
a uniform liquid is simply f3F ex = (n(p)) where V is the total volume. Santos introduces a correction so that 
$ = + &s the form of which is fixed by the scaling relations to be 

2s ( S 2 -3v 2 ) $ 3 (rj) \ 

(S 2 - V 2 ) $2 (f?) J 



$ s (n) = (s 2 - v 2 ) $ 2 [ri) A 



Here, the function A (y) is chosen so that the free energy in the uniform limit agrees with some chosen form (such 
as Carnahan-Starling 0). Calculation of the implied DCF in the uniform state based on Eq.(l) is straightforward 
(details are given in the Appendix). The result is 
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FIG. 1. The direct correlation function as calculated using the Percus-Yevik approximation, full line, the White-Bear approx- 
imation, dotted line, and Eq.Q with A(y) = (1 — (y/2) — log(l + y)/y)/3^. The left panel is for a low-density fluid with 
packing fraction r\ = 0.1 and the right panel is for a high-density fluid with r\ = 0.5. At low densities, the Percus-Yevik result 
is quite accurate and at high densities, the White-Bear result is in good agreement with simulations @|. 
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where x = r/cr and r\ — 5 per 3 and cpy (r; /?) is the Percus-Yevik DCF that comes from One feature that stands 
out is that this function diverges for x = unless A" (rz^j = 0- This divergence is unphysical and, as shown 

in the Figure, spoils the otherwise reasonable agreement with the results of the White-Bear functional. Hence, if 
this undesirable behavior is to be avoided, the only possibility is the relatively restricted set of corrections given by 
A (y) — a + by. However, Santos notes that in general one expects that A(y) ~ 0(y 2 ) so this eliminates the possibility 
of a correction. Santos also offers a modified version of his proposal that appears to avoid the divergence in the 
DCF but at the cost of violating his self-consistency condition]!]. It is interesting to note that similar terms arise in 
deriving the DCF from the Rosenfeld part of the free energy, but they cancel thus leaving the (finite) Percus-Yevik 
result (for a demonstration, see, e.g., the Appendix). 

In summary, the proposal of Santos produces a divergent DCF in the uniform liquid. Nevertheless, the reasoning 
behind it seems sound and it is to be hoped that there might still be a way to exploit it so as to eliminate some of the 
ambiguities of the usual extensions of FMT while retaining their advantages, one of which is an excellent description 
of the DCF for hard-spheres. 
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CALCULATION OF THE DIRECT CORRELATION FUNCTION 
1. The functional derivatives 



In general, 



Sp (n) ~5p(r 



-|$(n(r))dr (6) 



9 ™i /n(r)<Wl) 



and 



= J (n(r))«jj (r-rx)dr, = — $ 



5p(ri) (5/? (r 2 ) 
so that the bulk limit is 

lim 



y $ij (n (r)) Wl (r - n) ^ (r - r 2 ) dr, = $ (7) 



5 2 pF ex 



= $ij (n) y Wi (r - ri) tuj (r - r 2 ) dr (8) 



p(r)^p (5/9 (n) (5/9 (r 2 ) 

= (n) Wi * Wj 

where the last line introduces a compact notation for the convolution in which the spatial arguments are suppressed 
and where n = lim p ( r )^ p n (r). 

2. Useful formulae 

For later use, I note that the convolutions are given by 

(9) 



w n * w n 


= 2irO (a 


- r 12 ) 




-lf(x + 2) 


w s * w v 


= 2nO (a 


- n 2 ) 
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x) 


w s * w s 


= 27r9 (a 


- ri2) 
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Ax 




w Vi w Vz 


= 2irO (a 


- r 12 ) 




r) 



where x = r\ 2 ja. The weighted densities in the bulk limit are 

47T (0 \ 3 



^^1(2) p= r a (10) 

s(r) ->4tt(|)' 
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3. Rosenfeld functional 

It is straightforward to see that the Rosenfeld functional gives 
<5 2 } 

lim — ea - = 5$" (77) w T , * + 2$i (77) w s * w v (11) 

+ $2 iv) (s 2 - « 2 ) w »j * w -n + 2$2 (77) (w s * w s - io Vj * + 2$' 2 (77) to,, * (2sw s - 2viW Vi ) 
+ $3 (77) s (s 2 - 3w 2 ) to,, * w v + 2$3 (77) (s 2 - 3v 2 ) w n *w s + 2$ 3 (77) sw n * 2 (sw s - 3viW Vi ) 
+ 2$3 (77) w s * 2 (sw s - 3?;,^) + (f> 3 (77) s2 (w a * w s - 3w Vi * w Vi ) 



and since Vi = in the bulk this becomes 

A m x l \ T l \ = s$ i fa) Wj ) * Wl J + 2$ 1 fa) w * * w 'i 



Making use of 



+ $2 (77) s 2 w ?) * w,, + 2$ 2 (??) (w a * w s - 4$2 (»?) * w s 

+ $3 (77) s (s 2 - 3u 2 ) * w), + 2$ 3 (77) (s 2 - 3v 2 ) w, * w s + 4<J> 3 (77) s 2 w v * 

+ 6<J>3 (v) s ( w s *w s - w Vi * w Vi ) 



$! (77) = - J- In (1 - »?) $i (r?) = — ~o 77— — r <&'/ (r?) 



TTCT 2 (1 - T]) 7R7 2 (1 _ jj) 



27R7 (1-77) 2ttct (1 - r,) 2 2tr7 (1 - ^ 

^3(r ? ) = i-^^^(r7) = ^--^-3^$ 3 '(r7) " ' 



gives 



or 



247r(l-^) 2 av " 247r(l_^) J aw/ 24tt(1-^ 



hm —— — — — — r = - 7 (.x - 1) (a; + 2) (1 - a) 

2 

+ - 3 77 » (a: - l) 2 (a: + 2) - + 6 L (1 _ ^ 

fa-1) 3 »7-l (v-lf 

2 fa-1) fa-1) 

3 77 



Um J2 ^ } - 1 V {2v ± 1)2 3 ^ fa + 2 ) 2 X j ( 2? ? + !) 2 
p(p™pjp(n)<5p(r 2 ) 2 (I-77) 4 ^ 2 (l-r/) 4 (1 - t?) 4 

which is the usual PY expression. 

4. The new term 

Next, we need the additional terms introduced by Santos. First note that 

<5/3$s 

Xpfaj = 2 ( SWs ~ ^ W «i) $ 2 (77) A 
+ ( s 2 - w 2 ) <i> 2 (77)^ A 
+ (s 2 - v 2 ) $ 2 (»?) (A 5 «; s + A Vi w Vi + A n w n ) 

So, neglecting some non-contributing terms proportional to Vi, 



hm _ — - — ty~7 — \ = 2 ( Ws * w s — * w vi ) ^2 (t?) A + 4s (w s * w^) <I> 2 (77) A 



5 2 pAF^ ] 
P (r)^p 6p{r{)8p (r 2 ) 

+ 4s<f> 2 (77) * (A s w s *w s + ArfWs * w v ) + s 2 $ 2 (77) w v * w v A 

+ 2s 2 <I> 2 (77) (A s 77J n * 77J S + A^77J^ * Wr,) 

+ s 2 $ 2 (77) (A ss 7tj s * w s + 2A sri w s *w v + 2A v 2w Vi w Vi + A vv Wr, * w n ) 
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Using 



and the bulk limits 



gives 



A = A 2 



A, 



+ 



s(s 2 -3v 2 ) $ 3 fa)' 
(s 2 - v 2 ) $ 2 fa) 



s [s - 6V 



A' 2 



s s 



3v 2 ) $ 3 fa) 



(S 2 - V 2 ) $2 fa) 



$3 fa) 

(s 2 -v 2 ) \^ 2 (n) 



A' 2 



S S — 6V 



$3 fa) 



(s 2 -V 2 ) $ 2 fa) 



(S 2 -V 2 ) 2 $2 fa) J 



$3 fa A A , ( 2 *(s 2 -to 2 ) gafaA 



(s 2 - V 2 ) $ 2 fa) / 

(s 2 -3« 2 ) $ 3 



(s 2 - V 2 ) $ 2 fa) 



(s 2 -v 2 ) 2 $2 fa) y 

! V+3^)$3W\ / s(s 2 -3^ 2 ) $3 (g) ' 

85 ( S 2_ V 2)3 $ 2(7? )] I ( S 2_ u2) 



(18) 
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A s 
A, 

A„2 

A.« 



A 2s 



, *3 fa) 
'*2 fa) 

a' fV 3 ^ 



*2 fa) 



^ 2 fa)/ v *2 (77) y 



$2 fa), 

$3 fa)\ 



4$sfa)\ A , / $3 fa) 

S $2 fa) / V $ 2 fa) 



. jg fa) 
'*2 fa) 



A" 2s 



, $3 fa) 

$2 fa) 



(19) 



5 2 f3AF, 



plry^p 5p (ri) 5 p(r 2 ) 



2 (w s * w s - w Vi * w Vi )$2 fa) A + 4s (w s * w v ) $ 2 fa) A 

+ 4s$ 2 fa) (A s w s *w s + A v w s * w v ) + s 2 $ 2 ' fa) w v * w v A 
+ 2s 2 $ 2 fa) (A s w^ *w s + A v w v * w v ) 



s 2 $ 2 fa) 2 



, jg fa) V 

'*2 fa)y 



(20) 



A"w s * w s + 2A s?J w s * + A vv w n * w v + s 2 $2 fa) A' 



8 $3 fa) , 
s $ 2 fa) 



or, rearranging a little, 

Hm 5 2 (3AF^ 
P (r)^p <5p (ri)5p (r 2 ) 



(2 (w s * w s - to,,, * w Vi ) $2 fa) + 4s (io s * to,) $ 2 I 7 ?) + Mi) * w^s 2 $ 2 ' (77)) A 



(21) 



"'" *" ( $2 fa) ) * "" ) A ' 



+ s 2 $ 2 fa) 2 



, $3 fa) 
'$2 fa) 



w s * w s A" + 2A S)) w s * 10,, + A m w n * tu,, 



Next, using 



(s 2 -i; 2 ) 2 V^2 (T7)/ J \ {s 2 -v 2 ) <S> 2 {v) / 



$2(V)J (??) / \ $ 2(??) 



(s 2 -7J 2 ) \ < 1>2 (t?) / / \ {S 2 -V 2 ) <J> 2 (J?) 

+ / 2 ^( g 2 -3, 2 ) / 43(3) y\ A , / « (, 2 - 3« 2 ) $3 w ^ 



(s 2 - v 2 ) \$2(v)J I \ (s 2 - v 2 ) $2(1) 

--((li)T— (Si)"- 



results in 



nm 



P (r)^p 5p{r{)5p (r 2 ) 
= (2 (w; s * w s - io„ 4 * w Vi )$2 (v) + 4s (w s * w v ) $ 2 (??) + w v * w >i s2$ 2 A 



,*3 (V) 



+ 4s$ 2 (»?) — (tu s * w s - w^w^) + 4s ( j t. s «•,, + s" j ( ^ t [ j + - ( - r — - ) j «\. - ir n 

Inserting the convolutions and the bulk value of s gives 
llm — — r = - ri 2 )if (ri 2 /cr) 
if (x) = 27ra (??) + 677 (1 - x) $ 2 (ry) + ^r? 2 (a: - l) 2 (a; + 2) $ 2 ' (7?)^ A 

$ 3(t?)^ , c „ f $3 (»?) Y j , 3 2 /'/$3('))V . l/*3(ij)V'V. ...2 



+ 48 "*<"»^ + 6 H5^J <1 - I) + 5" 2 ^5^J + 3^J J(«-D J (- + «)JA' 

Finally, using the explicit forms for <3> 2 (t?) and $3 (77) and noting that 

$3 (t?) = 1 o 
$2 (77) 12 1 — 77 
2 „ $3 (t?) = 7? 
$2 (t?) 1-7? 

results in the form given in the main text. 
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